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Abstract 

This paper presents a model of a cryogenic refrigerator that integrates a reverse-Brayton lower temperature stage with a 2-piston 
Stirling upper temperature stage using a rectification system of check valves and buffer volumes. The numerical model extends the 
isothermal Schmidt analysis of the Stirling cycle by deriving the additional dimensionless governing equations that characterize the 
recuperative system. Numerical errors are quantified and the results are verified against analytical solutions in the appropriate limits. 
The model is used to explore the effect of the rectification system's characteristics on the overall cycle’s behavior. Finally, the model 
is used to optimize the hybrid system’s design by varying the swept volume ratio and phase angle in order to maximize the 
refrigeration per unit of heat transfer in the recuperator and regenerator. 
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1. Introduction 

This paper presents an ideal model of a cryogenic 
refrigeration cycle in which a recuperative lower stage is 
interfaced with a regenerative upper stage through a 
system of check valves and buffer volumes. The resulting 
cryocooler has a number of potential advantages. A low 
temperature recuperative stage avoids the dual loss 
mechanisms that plague regenerative cycles: thermal 
saturation and void volume pressurization. In addition 
to the anticipated higher efficiency, the configuration 
can be easily thermally integrated with refrigeration 
loads that are distributed in space or temperature. Me¬ 
chanical and electrical integration is also easy due to the 
continuous flow nature of the recuperative cycle. 

There are a number of applications that might utilize 
this generic concept of rectifying the AC flow within a 
regenerative system into a DC flow for a recuperative 
system. One example that is currently under develop- 
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ment is a pulse tube upper stage integrated with a tur¬ 
bomachine-based, reverse-Brayton lower stage, shown 
in Fig. la configured to cool an array of low tempera¬ 
ture superconducting electronic circuits. This configu¬ 
ration represents an extremely attractive combination of 
refrigeration technologies because the recuperative sys¬ 
tem can be activated by a relatively large pressure dif¬ 
ference from the linear compressor. This allows the use 
of simple hydrostatic gas bearings in the rotordynamic 
system for the turbine. Another possible configurations 
include a 2-piston Stirling cycle integrated with a Joule- 
Thomson low temperature stage, shown in Fig. lb. 
Without a recuperative stage, the rectification system 
alone may represent a convenient method of integrating 
a pulse-tube system with a distributed or isolated ther¬ 
mal load. 

The initial step in the development of these hybrid 
cryocoolers is to attain a solid understanding of the 
effect of the rectification system on the overall refriger¬ 
ation cycle. Although separately Stirling and reverse- 
Brayton devices have been thoroughly studied, the 
exercise of integrating an oscillating-flow cycle with a 
direct-flow cycle is new. This paper describes an extension 
of the isothermal Schmidt analysis of a 2-piston Stirling 
cycle by the addition of a rectification system (check 
valves and buffer volumes) and a recuperative stage 
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Nomenclature 

a area, m 2 

c P specific heat capacity, J/kgK 

m mass flow rate, kg/s 

m mass, kg 

p pressure, N/m 2 

q heat transfer rate, J/s 

q heat transfer per cycle, J 

r cv fluid resistance of check valve, N s/kg m 2 
R gas constant, J/kg K 

S reduced dead volume 

t time, s 

T temperature, K 

V volume, m 3 

w power, J/s 

w work transfer per cycle, J 

A dimensionless area 

M dimensionless mass flow rate 

M dimensionless mass 

Q dimensionless heat transfer rate 

0 dimensionless heat transfer per cycle 

R dimensionless fluid resistance 

W dimensionless power 

W dimensionless work transfer per cycle 
a phase angle, rad 

Greek symbols 

[> ratio of compression to expansion tempera¬ 

ture 

X ratio of expansion to cold load temperature 
<5 constant of solution 


</!> crank angle, rad 

y ratio of specific heat capacities 
i] efficiency 

k ratio of compression to expansion swept 
volume 

A ratio of buffer to expansion swept volumes 

n dimensionless pressure 

8 constant of solution 

r cycle period, s 

Subscripts 

b reverse-direction, check valve closed 

B buffer 

c, C compression space 

Carnot Carnot efficiency 

charge charge pressure 

cv check valve 

d, D dead volume 

e, E expansion space 

error numerical error 

f forward direction, check valve activated 

HP high-pressure buffer or check valve 

L cold load temperature 

LP low-pressure buffer or check valve 

max maximum 

rec recuperative 

rg regenerative 

St typical, closed Stirling cycle 

t turbine 

T total 


buffer 




Fig. 1. Potential hybrid regenerative/recuperative cryocooler configurations: (a) pulse tube/reverse-Brayton for LTS electronics; (b) 2-piston Stirling/ 
Joule-Thomson. 
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(heat exchanger and turbine). The number of governing 
equations is increased; reflecting the additional mass 
that is stored in the high- and low-pressure buffer vol¬ 
umes. The system of equations is solved numerically 
in order to determine the cyclic steady state pressure 
variations, which can be used to calculate various heat 
and work transfers. The model is presented in dimen¬ 
sionless form in order to identify the meaningful phys¬ 
ical groups for this hybrid system. In addition to the 
typical dimensionless parameters used to characterize a 
Stirling system (the swept volume ratio, phase angle, 
reduced dead volume, and compression-to-expansion 
temperature ratio), the hybrid system requires dimen¬ 
sionless groups that characterize the rectification system 
(the dimensionless normal and leakage check valve re¬ 
sistances, and the buffer volume ratio) and the recu¬ 
perative system (the dimensionless turbine size and the 
expansion-to-load temperature ratio). 

The use of an isothermal model even in the presence 
of varying pressure limits the model’s predictive value. 
However, the purpose of the model is to provide a 
qualitative assessment of this new type of cycle and 
develop an understanding the coupling between an os¬ 
cillatory and continuous system. The numerical model is 
developed in the first section of the paper. In the second 
section of the paper, some important performance 
metrics are defined and then derived in terms of the 
dimensionless solution. In the third section, the numer¬ 
ical model is used to explore the effect of the rectification 
system’s characteristics on the cycle’s behavior. In the 
fourth section, an optimization technique is described in 
which the phase angle and swept volume ratio are varied 
in order to maximize the lower stage refrigeration per 
unit of heat transfer in the recuperator and regenerator. 
The results of this optimization are shown to agree with 
prior studies of closed Stirling systems in the appropri¬ 
ate limit. The fifth and last section summarizes some of 
the important results of this work. 


2. Development of the model 

Fig. 2 illustrates a 2-piston Stirling refrigerator with 
isothermal compression and expansion volumes and an 
ideal regenerator. A first-order model of this device was 
initially developed by Schmidt [1] and subsequently re¬ 
stated and expanded in many textbooks on Stirling en¬ 
gines and refrigerators, for example Walker [2] or Urieli 
and Berchowitz [3]. This analysis is an appropriate 
starting point for an investigation of a hybrid configu¬ 
ration because the characteristics of the Schmidt model 
are well understood allowing the effect of the additional 
components related to the hybrid configuration to be 
readily identified. The first order analysis also provides a 
solid foundation for more detailed analyses that will 
follow. For example, the assumption of isothermal 



Fig. 2. 2-piston Stirling cryocooler. 

compression and expansion will be relaxed to attain a 
more accurate, adiabatic model similar to that used by 
Finkelstein et al. [4], The adaptation of the expansion 
process from a 2-piston mechanism to a piston-displacer 
or even a pulse tube configuration can also be achieved. 
Imperfect regeneration and recuperation can be in¬ 
cluded through the use of an ineffectiveness based on a 
total shuttled mass, as in Tailor and Narayankhedkar 
[5], or through a more sophisticated numerical solution 
of the heat equation within the regenerator matrix as it 
is subjected to the mass flow variations required by the 
refrigeration cycle. These first order results can be sub¬ 
sequently corrected for various parasitic losses such as 
pressure drop, shuttle heat loss, etc. which are post- 
calculated separately and applied to the net results, as 
described by Smith and co-workers (for example in [6]) 
and implemented by many others to achieve good 
agreement between analysis and data [7]. 

Fig. 3 illustrates the hybrid cooler that is formed by 
attaching a recuperative stage to the Stirling system 
through a rectification system that consists of two check 
valves, oriented in opposite directions, and two buffer 
volumes. The temperature-entropy diagram for the 
process undergone by the reverse-Brayton stage is il¬ 
lustrated in Fig. 4. Notice that the recuperative heat 
exchanger and turbine are both assumed to be ideal 
devices. This implies that the turbine inlet temperature is 
equal to the load heat exchanger exit temperature. The 
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compression piston 


compression space (T c , p, V c ) 
low-pressure check-valve, rh 


regenerator. 

expansion space ( T E , p, VJ 
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high-pressure buffer| 
( Te'Pb.hp• Vb) 


turbine 


turbine mass flow rate 


high-pressure check-valve, w 


low-pressure buffer 
( 7 © Pb.lpi V b ) 


recuperator 



refrigeration heat exchanger 




Fig. 3. Hybrid 2-piston Stirling/reverse-Brayton cryocooler. 


Pb.hp Pb.lp 



only loss in the reverse-Brayton stage, when considered 
alone, is related to heat transfer through a temperature 
gradient within the refrigeration load heat exchanger. 
The regenerative stage, as modeled, is reversible due to 
the assumed perfect components and isothermal heat 
exchange processes. Any additional losses in the system 
can therefore be attributed to the rectification system, in 
the form of pressure drop through the check valves. 

The analysis presented here will conserve much of the 
nomenclature used by Walker [2] in order to facilitate 
comparison of the results. The expansion volume (V e ) is 
assumed to vary sinusoidally according to 


Te(0) = y • [1 + COS(»] 


( 1 ) 


where V E is the swept volume in the expansion space and 
c/) is the crank angle 


The volume in the compression space ( V c ) is given by 

*c($) = ' [1 + cos(0 - a)] (3) 

where k is the ratio of the swept volume in the com¬ 
pression space to that in the expansion space and a is the 
phase angle. In addition, there is some unswept, or dead, 
volume (K D ) in the regenerative system. The total mass 
contained in the regenerative system ( m rg ) can be ex¬ 
pressed as the sum of the mass contained in these three 
volumes. 



kb 

7d 


P- Ve 


{/?• [1 +COS (</>)] 


2 R ■ T c 
-h K • [1 COS (0 — Of)] 2 • iS*]- 


( 4 ) 


where S is a dimensionless dead volume (sometimes 
called a reduced dead volume), and /? is the compression- 
to-expansion temperature ratio. 
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See 


Vd-T c 
V y Tv' 


(5) 


d/«B,LP 

At 


YYli fhcv ,LP 


(14) 


P = 


Tc 

T e 


( 6 ) 


Eq. (4) assumes that the temperatures of the expansion, 
compression, and dead volumes (7b, T c , and 7b, re¬ 
spectively) are constant, the instantaneous pressure in 
the system (p) is spatially uniform, and the working fluid 
is an ideal gas with a gas constant R. 

In a typical closed, Stirling system (denoted by the 
subscript St), the observation that the total mass re¬ 
mains constant allows Eq. (4) to be rearranged so that 
pressure can be expressed as a function of crank angle 


PstW 


_ 2 • R ■ T c ■ m Ig _ 

Vb, - {P • [1 + cos(0)] + k ■ [1 + cos ((/) — a)] + 2 • S} 

(7) 


Eq. (7) can be further simplified to [2] 

2 • R ■ T c ■ m Ig 

P%Aq>) Ve ■ {P + k + 2 ■ S) ■ [1 + 5 • cos(</> - 0)] 
where 0 and <5 are solution constants 


where w b .hp and «z b ,lp are the mass of gas contained in 
the high- and low-pressure buffer volumes, and m t is the 
mass flow rate through the turbine. The buffer volumes 
are assumed to be at a constant temperature (7b, the 
temperature of the expansion space) allowing Eqs. 13 
and 14 to be rewritten as 


d/?B,HP 

d(j) 


R-Te- t 

-(w cv ,hp -m t ), 

Vb • 2 • n 


(15) 


d/>B,LP 

dc/> 


R ■ Te • t 
V B ■ 2 ■ k 


■ (m, - Wcv.Lp) 


(16) 


where V B is the volume of each buffer. 

The mass flow rates in Eqs. (12), (15), and (16), must 
be expressed in terms of system pressures in order to 
solve the system of equations. A linear relationship be¬ 
tween mass flow through and pressure drop across each 
check valve is assumed 


P — P B,HP 

tn cv.HP — - 

r'cv.HP 


(17) 


0 = tan 1 


k ■ sin (a) 
p + k■ cos(a) 


(9) 


^cv.LP 


Pb.lp P 

^cv,LP 


(18) 


S = 


\Jft 2 + 2 • p ■ k ■ cos(a) + k 2 


2 S 


( 10 ) 


In the hybrid system, the mass of gas contained in the 
regenerative portion of the system is not constant and 
therefore the continuity equation for this volume must 
include mass flow to and from the recuperative stage 

dm, g _ . . 

i — ^cv,LP ^cv,HP (A 1J 


where 7M cv H p and tm CV iLP are the mass flow rates through 
the high- and low-pressure check valves, respectively. 
Note that either of these flows may be negative in the 
event that the check valves leak when exposed to an 
adverse pressure gradient. Substituting Eq. (4) into Eq. 
(11) leads to 


{p ■ [1 + cos(^>)] + k ■ [1 + cos(0 - a)] + 2 • S} ■ — 

+ p ■ {~p ■ sin ((f)) — k ■ sin (cf> — a)} 

_i-R'Tc ... N 

— i/ (^cv,LP ^cv,HP/ 

n • V E 


( 12 ) 


The pressures in the high- and low-pressure buffers 
(T’b.hp and /> b ,lp) continuously vary in response to mass 
flows from the regenerative system and through the 
turbine. Mass conservation for these buffer volumes is 
expressed by 


d«ZB,HP 

dr 


^cv,HP Mi-, 


(13) 


where r cv H p and r cv LP are the fluid resistance of the high- 
and low-pressure check valves. These resistances depend 
on the direction of the pressure difference 


F cv,HP — 


F cv,LP — 


^CV,f 

^cv,b 

P > T’b.hp 

T’b.hp 

(19) 

^cv,f 

F cv,b 

T’b.lp > P 

T’B.LP < P 

(20) 


where r cv j- is the forward resistance of the check valves 
and r cvb is the resistance to back-leakage. The flow 
characteristic of a turbine can be approximated as a 
nozzle over a wide range of operating conditions [8], 
allowing the turbine mass flow rate to be expressed as 


w, = 


T’b.hp 
R ■ T 


\ 


• a t 



/ ^ \ 


cp • 2 • 7 l • 

1 / i?B,HP \ 


\ T’B.LP ) 

_ 


( 21 ) 


where a, is an effective flow area that characterizes the 
turbine, 7b is the cold load temperature, c P is the specific 
heat capacity at constant pressure, and y is the ratio of 
specific heat capacities. 

Substituting the expressions for mass flow rates given 
by Eqs. (17), (18) and (21) into the continuity equations 
in the system, Eqs. (12), (15) and (16), it is possible to 
obtain a system of three coupled, ordinary differential 
equations for the pressures within the system 
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dp 

d(P 


P' 


§ • sin(0) + ^ • sin {<p - a) 


x-RTr 

n-V E 


% ■ [1 + cosO)] + H • [1 + cos {<p - 


[PB,LP—P) _ lP-lB,Hp) 
r cv,LP r cv,HP 


a)] + 2 • 


Vp-Tc 

v e t d 


( 22 ) 


d/>B 


d <p 


R ■ T E - t I {p — Pb.hp) 
2 ■ n ■ V B | r CV Hp 


PB,HP 

R-T, 


■ a t ■ 


\ 


2 • cp • 21 


1 - 


/?B,Hp 

Pb,lp 


1—7/7 


(23) 


d/>B,LP R Tp-t 

d (p 2 ■ n ■ V B 

Pb.hp 


„ rr, ‘ a t ■ 

RT l 
(Pb,lp ~p) 

cv.LP 


\ 


2 ■ Cp ■ /] 


L ' 


1 


Pb.hp 

T’B.LP 


1 — 7/7 


(24) 


Eqs. (22)-(24) suggest a set of dimensionless groups 
to characterize the hybrid system. The dimensionless 
pressure (77), check-valve resistance (7? cv ), turbine flow 
area ( A t ), buffer volume (A), and cold load temperature 
(x) are defined: 


77 = —— 


Rev = 


/^charge 
^cv ' 


A t = 


RT e -z 
a t ■ t 


V E 


sJy-R-Tp 




K.= 


Vb 

V e 

Te 


(25) 

(26) 

(27) 

(28) 
(29) 


where p cha rge is the charge pressure. 

The dimensionless form of Eqs. (22) (24) are 

d77 mW + K-sin^-a)] 

d <p {/8 • [1 -t- cos(<^>)] + k • [1 + cos (<p — a)] + 2 • 5} 


d77B,HP 

d<p 


(77 — 77b,hp) ■ 77b,hp • X 

2 • 71 • 7? cv hp • A 2 • 7t • A 


2 

(7- 1) 


( 77b,hp 

\ -ff B.LP 


1 — 7 / 7 ’ 


(30) 


( 31 ) 


d77B,LP 

d <p 


(77 b , lp — 77) A t ■ 77 b , hp • X 
2 • 71 • i^cv.LP ‘ 2 2 • 71 • A 


2 

(y- i) 


1 - 

/ 77b,hp \ 


\ 77b,lp / 


(32) 


These ordinary differential equations are stepped forward 
in time using a Runge-Kutta fourth-order technique 
starting from an initial charge pressure (77 = 77 b ,hp = 
77 b ,lp = 1) until a cyclic steady state is achieved. Cyclic 
steady-state is defined as the condition where the max¬ 
imum change in any of the dimensionless pressures is 
less than le- 4 between its value at the start and end of 
a cycle. The mathematical model is implemented in the 
EES—Engineering Equation Solver [9] software. 


3. Performance metrics 

The result of the model described in the previous 
section is the prediction of the cyclic steady-state pres¬ 
sure variation in the regenerative cycle and the two 
buffer volumes. In this section, we reduce these pressure 
variations into some important energy transfers that 
characterize the system. 

In order to compute the performance of the regener¬ 
ative stage, the heat transfer into the expansion space 
( q E ) must be computed. Provided both the regenerator 
and recuperator are perfect, this heat transfer is equiv¬ 
alent to the work transfer from the expansion space 
(t-Ve). 


p2-ti 

qE = w E = / p{(p) 

Jo 



(33) 


These and all subsequently considered work and heat 
transfers are made dimensionless (W and Q) through 
comparison with the product of the charge pressure and 
the expansion volume 


Qe=W e = 


H’e 


Pcharge * ^E 


= ~^f 0 n(<p)-sin((p)d(p 


(34) 


The power consumption of the refrigerator is related 
to the work transfer into the compression space (w c ). 
This is non-dimensionalized (Wc) and computed in a 
similar manner. 
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/'charge ' kf 
_ 1 
— 2 

The performance of the recuperative stage is related 
to the heat transfer into the refrigeration load heat ex¬ 
changer (< 7 L ). This heat transfer is equal to the work 
transfer from the turbine (w t ) which is, in turn, related 
to the mass flow rate through the recuperative stage. The 
instantaneous, dimensionless mass flow rate through the 
recuperative stage (M,) is 


f>2-n 


n{(f>) ■ sin (<p — a) dtp 


(35) 


m t • t • R ■ Tc 

/'charge * if 


— P ' X ' At ' 77b,HP 


2 

(y — i) 


1 - 

/ 77 ^hp\ 


\ 7 ?b,lp ) 


(36) 


The total dimensionless mass processed by the turbine 
per cycle (M t ) is obtained by integrating Eq. (36) 


M t = 


m t ■ R T C 

/'charge * If 


1 

2 • 71 



(37) 


The instantaneous power produced by the turbine (w t ) is 
given by 


w, = w, • c P • 7f ■ 


,, \ l ~y/y 

Pb,hp \ 

/'b.lp ) 


(38) 


Eq. (38) implies that the turbine is reversible and that 
the working fluid is a perfect gas (i.e. an ideal gas with a 
constant specific heat capacity). The turbine power may 
be transferred to a higher temperature as mechani¬ 
cal work transfer in a rotating shaft and dissipated as 
fluid power, as described by Swift [10], or transferred as 
electrical power and dissipated as resistive heating, 
as described by Nellis et al. [11]. The dimensionless 
instantaneous turbine power (W t ) is 


W t = 


W t ■ T 
kf ' /'charge 


Mi ■ y 

(y-!) -P-x 


( 77 b,hp 
\ 77 b,lp 


i —y/y 


(39) 


The total work produced by the turbine per cycle (W t ) is 
obtained by integrating Eq. (39): 


W t = 


w t 

kf ' /'charge 


1 

2 • 71 



(40) 


If the recuperator is perfect, then the turbine power is 
equal to the low temperature refrigeration load (Q L ) 


Ql = 


<7l 


kf * /'charge 


= Wt 


(41) 


The Carnot efficiency (>7 Ca mot) of a Stirling cycle 
modeled by an isothermal analysis is necessarily equal 
to one because no irreversible processes are included. 
However, some entropy generation is introduced as the 
reverse-Brayton cycle interacts with a constant temper¬ 
ature load, even in the limit of perfect recuperation and 
isentropic expansion. Also, the rectification system may 
introduce entropy generation related to the pressure 
drop across the check valves. The Carnot efficiency of 
the hybrid cryocooler is defined as the ratio of the power 
required by a reversible system supplying the same re¬ 
frigeration loads at the same temperatures to the actual 
power. 


7 Carnot 


&-•(/?•*-!) + &■(/?-1) 
Qc - Qe - Ql 


(42) 


Finally, the heat transfers in both heat exchangers, the 
regenerator and recuperator ( q Tg and q rec ), are of interest 
as these quantities determine the overall amount of heat 
transfer area required and, by extension, the ultimate 
size of the cryocooler. The total mass of gas that is 
shuttled through the regenerator during each cycle is 
equal to the maximum mass of gas contained in the 
compression space during a cycle. At any instant, the 
mass of gas in the compression space (m c ) can be cal¬ 
culated according to 

m c = 2 V( R ^ ■ [1 + cos {cp - a)] (43) 


This mass is made dimensionless (M c ) according to 

me ■ R • Tr n ■ k . , . .. ,. d , 

M c = - — = • [1 + cos {<p - a)] (44) 

/^charge "E ^ 


The heat transfer per cycle in the regenerator is 
related to shuttling fluid from the compression space 
to the expansion space: 

q rg = max(wc) • c P • (T c - 7f) (45) 


The regenerator heat transfer is made dimensionless 
(Qrg) according to 


Qrg 


qrg = Me ,max * y (/*-1) 
/'charge * kf (y 1) fi 


(46) 


The total amount of heat transferred in the recuperator 
during each cycle is obtained by integrating the instan¬ 
taneous rate of heat transfer 


?rec = / thf Cy ■ (7f — 7f) d t 
Jo 


(47) 


This recuperator heat transfer is made dimensionless 
(Qrec) according to 


Qrec 


qrec _ Mt-y-(x- 1 ) 
/'charge ' V E (? ~ 1) ‘ P ' X 


(48) 


If the check valves in the hybrid system are shut (i.e. 
if 7?cv,f = Rev, b —► oo) then the results predicted by the 
numerical model described above agree with those 
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predicted by the Schmidt analysis, verifying this aspect 
of the numerical model. The dimensionless input power 
predicted by the numerical model converges to the 
Schmidt solution to within 0.1% in this limit when more 
than 20 integration steps are used. The model has also 
been verified by comparing the instantaneous mass in 
the system to the charge mass; for more than 100 inte¬ 
gration steps this numerical error is found to be less than 
0.01%. At least 100 integration steps are used through¬ 
out the remainder of this paper. 

4. Effect of the rectification components 

Fig. 5 illustrates the pressure waves as a function of 
crank angle for the hybrid cycle as well as the pressure 
wave predicted by Eq. (49) for the Stirling cycle limit 
(i? cv ,f = f? cv ,b —► oo) . The most obvious effect of opening 
the check valves and connecting the recuperative stage is 
a reduction in the amplitude of the pressure wave ex¬ 
perienced by the regenerative stage. The crank angle at 
which the pressure reaches its maximum value is also 
slightly reduced. In the recuperative stage, the pressures 
oscillate more gently related to the cyclic draining and 
filling of the buffer volumes. The pressure difference 
between the high- and low-pressure buffers drives the 
turbine and produces refrigeration at the cold load 
temperature. The rectification system components are 
responsible for changing the oscillating pressure in the 
regenerative system (77) into the (more) steady pressures 
in each buffer volume (77 b ,hp and 77 B , L p). These com¬ 
ponents include the check valves and the buffer volumes 
and are characterized by the dimensionless groups 7? cv f , 
7? CVi b, and X. The effect of these dimensionless groups is 


straightforward; the optimal rectification system will 
have large buffer volumes compared to the flow per cycle 
required by the turbine (X/A x —> oo), check valves that 
offer no resistance to flow in the on-state relative to the 
turbine itself (7? eVif 4 t —> 0) and infinite resistance in the 
off-state (7? cv , b —> oo). Fig. 6 illustrates the pressure 
waves as the rectification system approaches this ideal 
limit, very large buffer volumes and check valves that 
approach perfect behavior. The regenerative system 
pressure oscillates between the pressure in the high- and 
low-pressure buffer volumes and these pressures are 
nearly constant. There is little pressure drop across the 
check valves when they are activated. The Carnot effi¬ 
ciency of this cycle is high because the ideal rectification 
system has not introduced any losses to the overall cycle. 
The turbine power produced by this cycle is constant— 
the infinite buffer volumes present it with a constant 
pressure ratio. 

Fig. 7 illustrates the situation where the check valve 
resistance remains negligible but now the buffer volume 
is finite. The buffer pressures are no longer constant and 
therefore the turbine sees a time varying driving pres¬ 
sure. While the Carnot efficiency of the cycle in this limit 
is still high, the fluctuating turbine power may be 
problematic for a real system; depending on the inertia 
of the rotating assembly the speed of the turbine may 
vary and in any event the turbine will be forced to op¬ 
erate away from its design point for a significant portion 
of each cycle. Fig. 8 illustrates the magnitude of the 
cyclic variation in the turbine power normalized against 
its average power as a function of the buffer volume for 
different values of check valve resistance. At high values 
of buffer volume the turbine power variation is small 
due to the stability provided by the large amount of fluid 



Fig. 5. Effect of rectification system on the pressure variations for R cv f = 0.2, R cv y —> oo, /J = 3.0, k = 2.0, S = 0.5, ct. = 1.0 rad, X = 1.25, y = 5/3, 
X= 3.0, A t = 0.25. 
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Fig. 6. Pressure waves in cycle for near perfect rectification system (R cv yA t approaching 0, R cv y approaching oo, X/A t approaching oo) for [S = 3.0, 
k = 2.0, S = 0.5, 0 ! = 1.0 rad, y = 5/3, x = 3.0, A, = 0.25. 



Fig. 7. Pressure waves in cycle with near perfect check valves (R^jAt approaching 0, S cvb —» oo) but finite buffer volume (A/A t = 5.0) for [S = 3.0, 
k = 2.0, S = 0.5, a = 1.0 rad, y = 5/3, / = 3.0, ^4, = 0.25. 


capacitance. As the check valve resistance is increased, 
the turbine tends to run in very short spurts twice 
per cycle that correspond to the pressure within the re¬ 
generative system hitting its maximum and minimum 
values. 

Fig. 9 illustrates the situation where the buffer vol¬ 
umes are infinite but now the check valve resistance is 
finite. In this limit, the pressure in the buffer volume 
remains constant but the pressure in the regenerative 
stage must significantly over- and under-shoot the buffer 
pressures in order to provide the pressure difference re¬ 
quired by the check valves. The net effect is to reduce the 
efficiency of the cycle due to this additional loss mech¬ 
anism. Fig. 10 illustrates the Carnot efficiency as a 


function of check valve resistance for various values 
of buffer volume. At a very low value of check valve 
resistance they introduce no losses and therefore the 
efficiency limits to a value dictated only by losses in 
the refrigeration heat exchanger. At very high values of 
check valve resistance, the efficiency approaches unity 
because there is no mass flow rate through the valves 
and therefore no entropy generation in the recuperative 
stage at all. In the intermediate region, the degradation 
in efficiency is dependent on the buffer volume; large 
buffer volumes tend to degrade the efficiency because 
the buffer pressures are not able to follow the pressure 
within the regenerative cycle and therefore the valve 
pressure drop is higher. 
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Fig. 8. Fluctuation in turbine power as a function of buffer volume at various check valve resistances for S cvb —> oo, /? = 3.0, k = 2.0, S = 0.5, 
a = 1.0 rad, y = 5/3, x — 3.0, A, = 0.25. 



Fig. 9. Pressure waves in cycle with infinite buffer volume (X/A t —> oo) but imperfect check valves (R cv /A, = 0.05, R c ,y —> oo) for ft = 3.0, k = 2.0, 
S = 0.5, a = 1.0 rad, y = 5/3, x = 3.0, A, = 0.25. 


From the discussion above it is clear that the rectifi¬ 
cation system should possess a dimensionless buffer 
volume to dimensionless turbine flow area ratio, X/A t , 
that is large (in excess of 5) in order to allow the turbine 
to operate near its design point at all times. Also, the 
product of the dimensionless check valve resistance and 
the dimensionless turbine flow area should be small 
(below 0.05) so that the rectification system does not 
introduce significant losses. 

5. Optimization of the hybrid cycle 

The performance of the hybrid cooler is a function of 
nine dimensionless groups. 


f?Li Qe, - *7Carnot! 6recj 2rg 

= 3(P,x,S,R cv y,R cv , h ,X, K,a,A t ) (49) 

The temperature ratios between the two stages ([I and /) 
are likely to be specified according to the requirements 
of a particular application. The dead volume ( S ), check 
valve resistances (7? cv ,r and R cv y), and buffer volume (1) 
are all important to performance, but they are not de¬ 
sign parameters in that there is a clear optimal limit for 
each of these groups. In a multi-stage refrigeration sys¬ 
tem, the ratio of the upper- to lower-stage refrigeration 
loads (Qe/Ql) is likely to be specified along with the 
temperature ratios. When all other parameters are fixed, 
this ratio can be adjusted directly by varying the di- 
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Check Valve Resistance, R cv f A { 

Fig. 10. Carnot efficiency as a function of check valve resistance at various values of buffer volume for R cw y —> oo, P = 3.0, k = 2.0, S = 0.5, a = 1.0 
rad, y = 5/3, % = 3.0, A t = 0.25. 


mensionless turbine area (A t ), which allows the turbine 
to accept a larger or smaller amount of flow. As the 
dimensionless turbine area is increased, more heat is 
lifted in the recuperative stage and less in the regenera¬ 
tive stage. 

This leaves two “free” parameters that can be varied 
to optimize the hybrid system: the swept volume ratio 
(k) and the phase angle (a). The isothermal Schmidt 
model of the closed Stirling cycle has been well studied 
and Walker [2] showed that there exists an optimal 
phase angle and swept volume ratio that maximizes the 
refrigeration power per unit mass of gas in the system. 
This metric is less meaningful in a hybrid system because 


a large part of the fluid inventory is contained in the 
buffer volumes and therefore does not participate in the 
refrigeration process. A more important parameter is 
the heat lifted from the recuperative stage per unit of 
heat transfer in the regenerator and recuperator. The 
total heat transfer in these heat exchange devices will be 
proportional to their physical size and therefore this 
optimization will provide the most refrigeration for a 
fixed amount of heat transfer area. 

Fig. 11 illustrates contours of refrigeration load per 
heat exchanger load as a function of the swept volume 
ratio and the phase angle and clearly shows that there is 
an optimal swept volume ratio and phase angle for the 



Fig. 11. Contours of constant recuperative stage refrigeration load per heat exchanger load for A/A t = 5.0, R cv yA, = 0.02, R cw y 
S = 0.5, y = 5/3, x = 3.0, Q E /Q L = 5.0. 


oo, P = 3.0, 
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hybrid system that makes the best use of the available 
heat transfer area. Walker [2] showed that the optimal 
phase angle for a typical closed Stirling system lies be¬ 
tween 0.55 and 0.65 rad over a wide range of conditions. 
In the hybrid system this optimal value increases to near 
0.85 rad, consistent with the earlier observation that one 
effect of the rectification system is to move the peak of 
the pressure wave to an earlier point in the cycle. Ac¬ 
cording to Walker [2], the optimal swept volume ratio 
for a closed Stirling system is proportional to the tem¬ 
perature ratio (p) and ranges from about k = 5 at 
p = 10 to k = 1 at p = 3. In the hybrid system the op¬ 
timal swept volume ratio is much higher, consistent with 
the fact that the compressor now serves the dual purpose 
of energizing the recuperative and regenerative refrig¬ 
eration stages. 

Fig. 12 illustrates the optimal swept volume ratio and 
phase angle as a function of the ratio of the regenerative 
to recuperative refrigeration load. At very high values of 
Qe/Ql, the system limits to a closed Stirling cycle and 
the optimization process approaches the appropriate 
optimal values of phase angle and swept volume ratio. 
As more refrigeration is provided by the recuperative 
stage the optimal swept volume ratio increases in order 
to provide additional, rectified flow and the optimal 
phase angle rises gradually. In the opposite limit, Q\]Q\. 
approaching zero, all of the refrigeration is supplied at 
the cold temperature and therefore the optimal swept 
volume ratio approaches infinity. In an analysis that 
accounts for recuperator and regenerator losses, a finite 
swept volume in the expansion space will be desirable 
even in this limit in order to intercept these heat ex¬ 
changer losses at the expansion rather than the cold load 
temperature. 


6. Summary 

A hybrid cryogenic refrigeration cycle is described in 
which an oscillating flow from an upper stage, regener¬ 
ative cycle is rectified and used to supply a continuous 
flow to a lower stage, recuperative cycle. The poten¬ 
tial advantages of this type of cycle include avoiding 
low temperature regenerative losses and alleviating 
mechanical, thermal, and electrical integration issues. 
There are numerous applications that might benefit 
from this generic concept and therefore this paper pre¬ 
sents a first order analysis of a hybrid system in which a 
2-piston Stirling cycle with a perfect regenerator and 
isothermal compression and expansion spaces is inte¬ 
grated with a reverse-Brayton stage with a perfect re¬ 
cuperator and a reversible turbine through a simple 
rectification system that consists of check valves and 
buffer volumes. The model is presented as an extension 
of the isothermal Schmidt analysis and the numerical 
results are verified against this limit. 

The model is used to explore the effects of the recti¬ 
fication system on the overall system performance and 
also to optimize the hybrid system for a given set of 
operating conditions. Some of the important results in¬ 
clude 

• The ideal rectification system has infinite buffer vol¬ 
umes and check valves that present no resistance to 
normal flow and allow no reverse leakage; in this lim¬ 
it the turbine produces a constant power and the rec¬ 
tification system introduces no losses. 

• When the buffer volumes are finite, the flow through 
the recuperative stage oscillates. This causes a cor¬ 
responding oscillation in the turbine power during a 



Fig. 12. Optimal values of swept volume ratio and phase angle as a function of 0 e/£?l for X/At = 5.0, i? CVi f4 t = 0.02, S cvb —> oo, f! = 3.0, S = 0.5, 
7 = 5/3, X = 3.0. 
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cycle. The magnitude of this power oscillation is 
quantified and presented as a function of the dimen¬ 
sionless buffer volume and check valve resistance. 

• When the check valve resistances are finite, the pres¬ 
sure drop through the valves causes a reduction in the 
cycle efficiency. The magnitude of this loss is quanti¬ 
fied and presented as a function of the dimensionless 
check valve resistance and buffer volume. 

• There exists an optimal phase angle and swept vol¬ 
ume ratio that maximizes the heat lifted from the re¬ 
cuperative stage per load on the heat exchangers 
(regenerator and recuperator). The optimal value of 
the phase angle is slightly higher than in a “pure” 
Stirling system and the optimal value of swept vol¬ 
ume ratio is much larger. These values are presented 
for one particular set of operating conditions as the 
ratio of the heat lifted by the recuperative stage to 
the regenerative stage is varied. 
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